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The  Convergence  Of  The 
State  Probabilities  In  A  Class  Of 
m-Dimensional  Simple  Epidemic  Models 

by 

Herbert  Lacayo  and  Naftali  A.  Langberg 
ABSTRACT 

population  of  susceptible  individuals  exposed  to  m  contagious  diseases 
is  considered.  The  progress  of  this  epidemic  among  individuals  is  modeled  by 
an  m-dimensional  stochastic  progress.  The  components  of  this  process  represent 
the  number  of  infective  individuals  with  the  respective  diseases  at  tins  t. 

A  class  of  m-dimensional  stochastic  processes  is  constructed.  These 
processes  describe  the  progress -of  the  epidemic  models  considered  in  the 
sequel.  Exact  and  approximate  formulas  for  the  joint  and  marginal  state  prob¬ 
abilities  of  these  models  are  obtained.  It  is  shown  that  the  approximate 
formulas  are  very  simple  functions  of  time  while,  the  derivation  of  the  exact 
formulas  involve  tedious  computations.^. 

Key-words:  m-dimensional  simple  epidemics  and  stochastic  processes, 
exponential  and  negative  binomial  random  variables,  convergence  in  distribution. 
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1.  Introduction  and  Suag&ry. 


In  a  simple  epidemic  situation  we  assume  that  a  population  of  susceptible 
individuals  (suscootiblcs)  is  exposed  only  to  one  contagious  disease  (disease) 

C Bailey  (1975)3.  However,  frequently  susceptibles  are  exposed  simultaneously 
to  more  than  one  disease,  as  is  the  case  with  different  types  of  flu.  In  this 
paper  a  population  of  susceptibles  exposed  to  m  diseases,  m  «*  1,  2,  ...,  is 
considered.  We  say  that  this  population  undergoes  an  m-diaensional  simple 
epidemic  if  the  following  six  assumptions  hold. 

(1.1)  At  each  point  in  tine  at  most  one  susceptible  contracts  a  disease. 

(1.2)  Each  susceptible  contracts  at  most  one  disease. 

(1.3)  Once  a  susceptible  contracts  disease  r,  r  =  1,  . ..,  m,  he  remains  conta¬ 
gious  for  the  duration  of  the  epidemic. 

(1.4)  An  infective  individual  (infective)  with  disease  r  r  ®  1,  ...,  m,  can 
transmit  only  that  disease. 

(1.5)  All  interactions  between  a  susceptible  and  an  infective  with  a  specified 
disease  are  equally  likely  to  result  in  an  infection. 

(1.6)  Individuals  neither  join  nor  do  they  depart  from  the  population. 

For  m  *  1  the  m-dimensional  simple  epidemic  models  reduce  to  the  traditional 
single  epidemics. 

Let  Tq  denote  the  first  time  we  have  at  least  one  infective  with  each  of 

the  various  diseases,  and  n  tho  number  of  susseptibles  at  time  T^.  The  progress 

of  a  m-dimensional  simple  epidemic  among  susceptibles  is  described  by  an  m-dimen- 

sional  stochastic  process  X  (t)  -  [x  (t),...,X  (t)],  teCO,-).  The  components 

“*i»  n  f  x  n  |Ni 

of  X^(t)  represent  the  number  of  infectives  with  the  respective  diseases  at 
time  t  measured  from  TQ.  In  Section  2  a  class  of  m-dimensional  stochastic 


2 


processes  is  constructed.  These  processes  describe  the  progress  of  the 
epidemic  models  considered  in  the  sequel. 

The  computation  of  state  probabilities  in  epidemic  models  is  of  major 
interest  to  researchers.  In  Section  3  we  derive  formulas  for  the  joint  and 
marginal  state  probabilities:  P_  .  (t)  *  P(X  (t)-X  (0)=k},  and  P  .  (t)  * 

U|K  “'ll  il  “* 

P{Xn  r(t)-Xn  r(0)=k>,  where  k  »  [kj,...,^],  k,  kj,  ...»  kme{  0,1,...}, 
r  ■  1,  ...,  m,  and  te[0,°») .  These  formulas  are  calculated  without  the  tradi¬ 
tional  use  of  the  differential  equations  associated  with  the  state  probabilities. 
This  is  done  by  utilizing  a  formula  for  the  distribution  function  (df)of  a  sum 
of  independent  exponential  random  variables  (rva's)  given  by  Billard,  Lacayo 
and  Langberg  (BLL)  (1980) . 

-  -  *  ■  •  •  •  m 

The  formulas  for  the  state  probabilities  obtained  in  Section  3  are 
rather  complicated .  To  overcome  this  deficiency  we  derive  very  simple 
approximations  to  the  joint  and  marginal  state  probabilities  when  the  initial 


number  of  susceptibles:  n,  is  sufficiently  large.  Let  us  denote  X  (0)  by 

Hir 

b  ,  r  ■  1,  ...»  m.  In  Section  4  it  is  proven  that  for  all  t€(0,-),  and 


r  ■  1, 
(1.7) 


f"  \k(t) 


m  fb+k  -1  -o  b  t  -at  k 
n  l  l  e  r  r  (1-e  r  )  r 
r-1  br_1  . 


(1.8) 

(1.9) 


b  *k-l]  -a  b  t  -a  t  k 
Pn  k  r(t)  "  b  -1  ®  (1'e  >  * 

ft*-  n,K»  l  r  J 

at 

tin  EX  ft)  =  b  (e  -1),  and  that 

ft*-  n,r  r 

2b  t  at 

am  Var(X  ft)}  =  b  fe  r  -e  r  ). 

n y &  r 


(1.10) 


BLL  (1979)  consider  a  special  class  of  m-diraensional  simple  epidemic 
models  and  name  them  the  symmetric  m-dimensional  simple  epidemics.  These 
models  are  a  subclass  of  the  ones  considered  by  us,  and  are  presented  in 


Section  2.  BLL  (1979)  prove  Statements  (1.7)  through  (1.10)  for  the  symmetric 
models.  The  results  proven  in  Section  4  generalize  those  obtained  by  BLL  (1979) 
Further,  since  the  methods  used  by  BLL  do  not  apply  to  the  more  general  case 
considered  by  us  a  new  way  of  attacking  the  problem  is  presented. 

Finally  we  mention  that  the  process  X  (t)  can  be  used  to  describe  various 

*il 

competative  situations  other  than  the  m-dimensional  simple  epidemic.  Two  such 
competative  situations  are  presented. 

(I)  In  marketing,  let  X  (0)  represent  the  numoer  of  customers  who 

n,r 

initially  own  brand  r  of  a  product.  Let  n  be  the  number  of  customers  who  at 

the  time  we  start  to  observe  do  not  possess  the  product.  Then  X  (t)  can 

n*r 

represent  the  random  number  of  customers  who  own  brand  r  at  time  t. 

(II)  Consider  an  election  campaign.  Let  X  (0)  denote  the  number  of 

nfr 

electors  who  support  party  (or  candidate)  r  at  the  beginning  of  the  campaign; 

and  let  n  be  the  number  of  unconnitted  voters  at  that  time.  Then  X  (t)  can 

nfr 

describe  the  random  number  of  voters  that  support  party  (or  candidate)  r  at 
time  t. 

The  various  brand  manufacturers  (political  parties  or  candidates)  are  competing 
for  customers  (voters)  as  diseases  are  "competing"  for  susceptibles .  Thus, 
the  results  obtained  in  this  paper  are  of  interest  not  only  to  researchers  in 
Epidemiology  but  as  well  to  those  who  deal  with  various  competitive  real  life 
situations . 


n--.  1r  r  •'  rlil»tii%rt 


2.  Model  Construction. 

In  this  section  we  construct  a  class  of  m-dimensional  stochastic 

processes.  These  processes  describe  the  progress  of  m-dimensional  simple 

epidemic  models  considered  in  the  sequel. 

Some  notation-  is  needed.  Let  X  ft)  »  J®  .  X  ft) ,  be  the  total  number 

nk  _  tr=l  n,rv 

of  infectives  at  time  t  measured  from  Tp ,  te[0,~) ,  and  Tr  ^  the  k—  inter¬ 
infection  time  defined  as  the  time  that  elapses  between  the  (0)+k-l  and  the 
X  (0)+k  infection,  k  ■  1,  ...»  n.  Further,  let  £  .be  a  rva  assuming  values 
in  the  set{ 1 , . . . ,m}  designating  the  disease  responsible  for  the  Xn(0)+k  in- 

_Jr 

fection,  k  »  1,  ...»  n,  and  let  £  _  =  0.  Finally,  let  S  .  =  J  _»T  „  be 

n |U  R|K  (jsi  n 

the  time  measured  from  Tp  until  the  X^ (0)+k  infection,  k  =  1,  . ..,  n,  Sn  p  =  0, 

S  _  -  ■  •»,  and  let  I  be  the  indicator  function.  In  particular  S_  ^  is  the 
n,n+i  —  n,n 

duration  time  of  the  m-dimensional  single  epidemic. 


For  k  =  0,  ...»  n,  r  *  1,  ...,  m,  and  te(0,»),  the  following  event  equality 


holds . 


(2.1) 


lXn,r<«-Xn.r»»*k>  = 


Thus,  to  construct  the  process  X^(t)  it  suffices  to  determine  the  df  of  the 

random  vector  (rve)  [T  ,£  ,q»l,...,n).  Before  determining  the  df  of  this 

-  n,q  n,q 

rve  we  introduce  some  more  notation.  Let  J  ,  _  be  the  index  set  of  all  in- 

n»k»r 

fectives  with  disease  r  at  time  TQ+Sn  and  let  cn  ^  r  *  ^n  q“r) 

be  the  number  of  infectives  with  disease  r  at  time  Tp+Sn  ,  k  «  1, . . . ,n, 

r  *  1,  ...,  m.  Further,  let  x  f.  .  bo  a  rva  that  describes  the  time  measured 

th 

from  the  Xn (0)+k-l  infection  until  the  i—  contagious  individual  causes  the 
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Xn(0)+k  infection,  i  -  1,  Xn(0)+k-l,  k  *  1,  ....  n,  and  let  <^,...,0^(0,-), 

Throughout  we  assume  that 


(2.2)  H>.  conditional  r*a's  hn#i>klt,<0 . (n  k.j),  i  -  1 . X„CO).k-l. 

k  *  1,  n,  are  independent  exponentially  distributed,  and  that 

(2.3)  E(tn  ±  k|5n >0»**.#Cn ,k-l)*nC«rO»-k*l)r1  for  k  =  1,  ...,  n, 

r  a  1,  •••,  n. 


We  are  ready  to  determine  the  df  of  the  rve  CT  ,K  ,q=l,...,n).  The 

n,q 


following  two  lemmas  are  needed. 

Lemma  2.1.  Assume  Conditions  (2.2),  (2.3)  hold.  Then  for  r  =  1,  ....  m, 
and  k  *  1,  ...,  n 


(2’4)  P{€n,k”r^n,0’*“»^n,k-1}  " 

*  °rC n,k,r^»l°£Cn,k,£^  * 


Proof.  Let  k,  r,  be  fixed.  Then  the  event  {£  ,»r)  is  equal  to 

nvK  ^ 


m 


("1"CTn.i.k!ieJn,k,r]<"d"^n,i,k:1'“iiIJ„.lt,z2>- 


**r 


Consequently  Statement  (2.4)  follows  by  Conditions  (2.2),  (2.3)  and  some 
simple  integral  evaluations .  | | 

Let  «(l,...,m),  and  lQ  »  0.  Then 


(2.5) 


^n.c-V-1 . n)  ' 


n 


'  . ,-U- 


Thus,  Statement  (2.4)  determines  the  df  of  the  rve  [£  C  3. 

n,n 

Lemma  2.2.  Assume  Conditions  (2.2),  (2.3)  hold.  Then  the  conditional 


\ 


•  ,  ■  ... 
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"a'S  «Vj*M . W>-  q  ■  1,  ...,  n,  are  independent  exponentially 

distributed  with  means  respectively  equal  to  [(n-q+lHE^a^  q  r)n_1  ]_1 . 

Proof.  Note  that  Tn  ^  min  [tn  .  q:  i-l,...,Xn(0)+k).  q  «  1,  ...,  n. 
Consequently  the  result  of  the  lemma  follows  by  Conditions  (2.2),  (2.3),  and 
some  single  properties  of  exponential  rva's.  || 

Clearly  Lemmas  2.1,  2.2  together  determine  the  df  of  the  rve 
*  *  *  *n^* 

Although  the  transition  rates  of  the  various  diseases  are  not  used 

explicitly  they  are  presented  for  the  sake  of  completeness. 

Definition  2.3.  The  transition  rate  of  disease  r  at  time  t,  r  »  1,  ...,m, 
t€Co,»»),  is  given  by  ....  . 

-tin  h_1P(X  (t+h)-X  (t)=l|X  (t)}, 

fc*0+  ’  n,r  ^ 

and  is  denoted  by  R(Xn(t),r). 

Finally  it  is  shown  that  the  transition  rates  of  an  m-dimensional  simple 
epidemic  that  satisfies  Condition  (2.2), as  expected,  determine  uniquely  the 
epidemic.  We  need  the  following  lemma. 

Lemma  2.4.  Assume  Conditions  (2.2),  (2.3)  hold.  Then  for  r  ■  1,  ...,  m. 


and  te[0,») 


(2.6) 


R(X„(t),r)  -  n‘*ar  Xntr(t)(n“Xn(t)+Xn(0))* 


Proof.  The  results  of  the  lemma  clearly  follow  by  the  memoryless  property 
of  exponential  rva's  [Barlow,  Proschan  (1975),  p.  56),  Equation  (2.1)  and 
Conditions  (2.2),  (2.3).  || 

Consequently  we  obtain  for  ieJ  .  ,  k  ■  1,  ...,  n,  and  r  »  1,  ...,  m,  that 

JifK  fr 


(2.7) 


EtTn.i,kl5n,0#**,,5n,k-lJ  * 

"  Cn,k,r{R(CCn,k,l*’-*»Ci».k,in:,»r))  *• 


Thus,  m-dimensional  single  epidemic  models  that  satisfy  Condition  (2.2)  are 
uniquely  determined  by  their  transition  rates. 

BLL(1979)  define  a  class  of  m-dimensional  simple  epidemic  models  and 
name  them  the  symmetric  m-dimensional  simple  epidemics.  Specifically  a 
population  of  susceptibles  exposed  to  m  diseases  undergoes  a  symmetric 
m-dimensional  simple  epidemic  if  Assumptions  (1.1)  through  (1.6)  hold.  Condi¬ 
tion  (2.2)  is  satisfied,  and  the  transition  rate  of  disease  r  at  time  t, 
r  *  1,  ...»  m,  te(0,«),  is  given  by 

(2.8)  RQ^CO.r)  »  r(t)(n-Xn(t)+Xn(0)),  where  <*€(0,~). 

The  transition  rates  given  by  Equation  (2.6)  reduce  to  those  of  Equation  (2.8) 
by  selecting  ...  *  am  «  o.  Thus,  the  symmetric  m-dimensional  simple 

epidemic  models  are  a  particular  case  of  the  ones  constructed  in  this  section. 
For  the  symmetric  m-dimensional  simple  epidemic  we  obtain  from  Lemma  2.2  that 
the  interinfection  times:  Tn  j,  . ..,  Tn  n,  are  independent  exponential  rva's. 
independent  of  the  infection  causes:  Cn  •••»  Sn  n*  This  special  structure 
does  not  hold  in  general  for  m-dimensional  simple  epidemic  models.  BLL  (1979) 
use  this  special  structure  to  prove  Statements  (1.7)  through  (1.10)  for  the 
case  of  symmetric  m-dimensional  simple  epidemic  models.  Their  method  of  proof 
does  not  apply  to  the  models  considered  by  us.  To  prove  Statements  (1.7) 
through  (1.10),  for  the  models  constructed  in  this  section,  a  new  way  of 
attacking  the  problem  is  needed.  One  such  way  is  presented  in  Section  4. 


^>JS,  -  A-  -  —  *-■»- 
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3.  Formulas  for  the  State  Probabilities. 


This  section  contains  formulas  for  the  joint  and  marginal  state 
probabilities:  P  v(t)  and  P  .  (t),  r  =  1,  . ..,  m,  te(0,«).  These 

formulas  are  calculated  without  the  traditional  use  of  the  differential 
equations  associated  with  the  state  probabilities.  Rather,  we  utilize  an 
available  formula  for  the  distribution  function  of  a  sum  of  independent  ex¬ 
ponential  rva's.  For  the  sake  of  completeness  this  formula  is  presented. 

We  need  some  notation.  Let  0^,  &2>  •«*»  he  positive  real  numbers,  and 

Uj ,  U2,  ...,  be  iid  exponential  rva's  with  mean  equal  to  one. 

Theorem  3.1.  CBLL(1930),  Theorem  23.  Let  M  be  a  positive  integer, 

A,, Cl)  •  t-l)M-‘  ’  9  and  *^(j)  -  (-1)H  l  “  s„q  .  j  -  M  .  1, 

^  q-i  “  *  1-1  ’ 

M  +  2 .  Then 


(3.1) 


P{JM  .e-1u  St) 
^qal  q  q 

r 


*  2j*M(j!)"  VJ)(-tr»  ««©*-). 

To  aid  in  computing  the  joint  and  marginal  state  probabilities  some  more 
notation  is  introduced .  Let  L^t  . . . ,  be  a  sequence  of  positive  integers 
assuming  values  in  the  set  {1 . m),  lQ  •  0,  and  let  s(k)  »  .  Further, 

let  B(k)  «  {C^,...,£s(yJ:  Iq.(^U^q»r)“kr,r»l. ••.,»),  and  Afr.j.k)  - 
{ , . . . 3 :  (£q*r)®k)  for  j  •  k,  k  t  l,  .»»,  r  *  1,  ,  m,  and 

k  -  0,  1,  ...  .  By  Equation  (2.1)  we  obtain  for  te(0,«)  and  r  -  1,  ...,  m,  that 


(3.2) 


u\ . .  . *<*» 


QS) 

#  .Ms(y )  I 


* 

'■'V! 

£ 


*s(y)«B<!c)). 
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<3-3>  pn,l..rC*>  - 

•  ln  l  tPtSn,jst<Sn,j+l^n,qs'^q»q*0,***^> 

j*k  Uj 

*  PUn ,qa£q »qs*0 » *  *  * » } 1  tt£i .  *  •  •  3eA (r , j  ,k) )  ] . 

Thus,  to  compute  Pn  k(t)  and  Pn  k  r(t)  it  suffices  to  evaluate  PUn^®£q»q*0» 

and  P{Sn  kst<sn>k+lUn,q**q»t»3°»‘--»kJ  for  k  “  °»  •••»  n*  Let  »0  "  l»  and 
u  *  ?m  a  [b  ♦Fq”*I(£.»r)],  q  *  1,  2 .  For  simplicity  in  notation  the 

dependence  of  yq  on  £ql,  q  *  1,  2,  ....  is  suppressed.  First 

P{E  =1  ,q=0,...,k}  is  evaluated. 
n,q  q'M 

•  «  •  •  ■  Lemma  3.2.  .Let  ke{.0,l,.«},a  and  d^.  *  Iq3Q^^qs^)  • 


Then 

(3.4) 


P{?n,q=£q^° . 10  * 


~  m  d 

fb  +d  -ill 

n  d  !a  r 
.  r  r 

r  r 
b-1 

r=l 

r  Ji 

k  -1 
(  I-1]. 

q=0  H 
k 

Proof.  Note  that  P{£n  q®£q»q®0» . . .  ,k}  =  II  P(£n^q®-£q|£n^  ja^j  » 3S®»  •  •  •  • 

Consequently  Statement  (3.4)  follows  by  Equation  (2.4).  || 

Next  P(Sn  k^t<Sn  ^1|£nq»iq,q®0,...,k>  is  computed.  Let  fq(61> . . . ,0q.t)  be 

the  density  function  of  the  rva  JqxjO^Uj,  q  *  1,  2,  ...  .  Then  for  k  **  0,  1,  ..., 

and  te  (0  ,*•)  . 


Thus ,  by  Lemma  2 . 2  and  Theorem  3 . 1  for  k  *  0 ,  1 ,  . . . , 


and  te (0,®) 


(3'6)  HS„,kit<S„,kJ5n,q*V‘>‘° . 10  * 


•Ujt 


k  =  0 


r  uii-Vf-i)"*1 


n 


j*n+l 


j1+--+jn=j  q=i 


l  n  |Tq  +(-l)n_1  (n!)_1tn  n  U  k  «  n 


.1  «l 


q=i  q 


k+1  j-k+2  J1  •••  3k*l+j  q*l  q 


k+l  j  +1 


Uk 


-1  k  -i^+* 

*  »*!/(*!)  n  » 

q=l 


The  formulas  for  the  joint  and  marginal  state  probabilities  are  obtained  from 
Equations  (3.2)  and  (3.3)  by  subtitution.  Note  that,  as  expected,  the.  exact 
formulas  for  the  joint  and  marginal  state  probabilities  are  rather  complicated 


-*?>:••  -  -  ,  ■  . 

„.«*«■  -•  VW-,-  - 
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4.  Asymptotic  Approximations . 

In  this  section  we  present  the  asymptotic  approximations  of  the  joint  and 
marginal  state  probabilities  and  some  related  moments.  All  limits  are  calcu¬ 
lated  as  n  ■*  «. 

First  it  is  shown  that  the  joint  and  marginal  state  probabilities  satisfy 
Statements  (1.7)  and  (1.8).  Three  lemmas  are  needed.  Throughout  these  lemmas 

let  \  a  ^r»lor^j-0I(£j“r)»  and  **  in  Section  3  let  p0  =  1,  and  uq  *  V^=l°\br' 
q  =  1 ,  2,  ...  .  For  simplicity  we  suppress  the  dependence  of  and  on  t.b* 
sequence  ^  $  ...  . 

1. 

Lemma  4.1.  Let  ke(0,l ,...},  and  dr  =  Iq-Q1^*1) »  r  *  1» 

(4.1) 


m.  Then 


£im  P( n-l  ,q=0 , . . .  ,k}  = 
n,q  q 


l 


a 


*=  n  d  !a 
r=l  r  r 


b  +d-l 
r  r 

b  -1 
r 


C  n  u  ] 

q*0  q 


-1 


J  <1 


Proof.  The  desired  result  follows  clearly  from  Lemma  3.2.  (| 

Lemma  4.2.  Let  ke(l,2,...}.  Then  the  conditional  rva's 

Zk  —1 

,W  u  . 

«,H  M  q»l  q  q 

Proof.  To  prove  the  desired  result  it  suffices,  by  the  Cramer-Wold  device 

[Billingsley  (1968),  p.  49],  to  show  that  the  conditional  rve's 

{[T  -,..., T  .  j|C  *£  ,q=0 .... ,k-l)  converge  in  distribution  to  the  rve 

Hjl  n,K  n*q  q 

[p”^Uj, . . . .v"1!^].  The  preceding  statement  follows  by  Lemma  2.2.  || 

To  prove  the  result  of  the  next  lemma  we  need  [Renyi  (1970),  p.  203], 
the  following 


gflWVl  ■  fN(*l . V‘>  •  iJJ.lW 


■X  t 

q 


(4.2) 
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N 


,-l 


where  Aj,...,X„.  aro  distinct  positive  real  numbers  and  C  (q)  -  [  n  X.][  ji(A  -X  )]** 

-  N  j-1  3  j»l  3  q 

i  3+n 

q  »  1,  ... t  N. 

Lecuna  4.3.  Let  kj,  ....  kffl  «(0,1,...},  k  »  Ck|,...,k  3,  and  H(t)  = 

m  k  txs  i  -V.t  s(y+i 

(  "  k  !®  >  l  I?=r  ®  3  C  n  Cv.-v.jj”1.  teCO,-).  Then 

r=l  r  r  UfV\  j  1  i-i  1  J 


(4.3) 


B(k) 


i=l 

m  -at  k 

H(t)  =  n  (l-e  r  )  ,  te[0,«). 

r*l 


Proof.  The  proof  is  by  an  induction  argument  on  k. 

First  note  that  Statement  (4.3)  holds  for  k  =  1 .  Assume  Statement  (4.3)  holds 

for  k  -  1  (k>J).  Let  rq  -  JJ.2  2J.1«rI(V)'  q  ’  2'  5 . •«*  • 

. V-  C-£j .... jeB(k) } .  Now  by  the  induction  assumption 

(4.4)  dH(t)  , 
dt 

_  -at  k  -1  m  k.  ,  .  -y.t  s(k)+l  . 

“  5-1®  kr0rC(V1)!or  n  ki°i  l  S«2e  C  n  <VV3  * 

T  1  *  r  r  r  i*l  1  1  B(r,k)  3 z  i-2  1  3 

i#r  i+j 


-at  k 


-a  t  k 


_  -at  .  at  m  -at  k  m  -»«.*. 

*  5-lkr°r(1"e  )  le  n  (1_e  >  T  ■  4J  n  (J-e  r  )  T> 

r-1  dt  r»l 

By  Equation  (4.2)  Ptl^jX^U^st)  -  l-^jelX"1CN(q)e  q  .  Since  Uj,  v2»  are 

distinct,  we  conclude  by  (4.2)  that  H(0)  »  0.  Consequently  Statement  (4.3)  for 

k  follows  by  integrating  both  sides  of  Equation  (4.4).  || 

We  are  now  ready  to  prove  Statement  (1.7). 

Theorem  4.4.  Let  t«(0,»).  Then 

m 


*i®  Pn  k(t)  -  n 
n»K  r*l 


W1 

V1 


e  “^(l-e"*^)1*. 
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(4.5) 


Proof.  By  Equation  (3.2)  and  Lemmas  4.1,  4.2  we  obtain  that 
)  p„  k(t)  -  l  CP(L  Vn'1U  it<^f^+iu"1u  }  • 

tfl—*  1  q  q 

f  "  k-fbr+I;r-1]]  s(]i)  1 

*i  Var  !  J  -t  n  ([£.,..., ^JeBCk))] 

(r-l  l  r  -  JJ  q«l  **  1 


By  Equations  (3.5)  and  (4.2) 
(4.6)  £im  Pn  k(t) 


Consequently  the  desired  result  follows  by  Lemma  4.3.  J  j 

For  the  sake  of  completeness  we  present  the  following 

Definition  4.5.  We  say  that  W  is  a  negative  binomial  rva  with  parameters 

ae{),2,...}  and  pe(0,l),  and  write  W^NB(a,p)  if  P(W»k}  »|a**j1jpa(l-p)*,  k«0,l,.. 

Next  Statement  (1.8)  is  proven.  Let  X_(t)  ■  CXj(t),...,Xm(t)),  t«(0,«*)  be  a 

collection  of  rve's  with  independent  components  such  that  Xr(0)  *  bf,  and 
'art 

Xr(t)M<B(br,e  ),  r  =  1,  te(0,»).  Since  ^(t)  and  X(t)  are  discrete 

rve's  we  obtain  by  Theorem  4.4  and  a  well  known  result  (Billingsley  (1968),  p.  16 


that 


Cor'*1'* «trv  4.6.  Let  tc (0,») .  Then  the  rve's  X_(t)  converge  in  distribution 
to  the  rve  £(t). 

By  the  Cramer-Wold  device  and  Corollary  4.6  we  conclude 

Corollary  4.7.  Let  te(0,»)  and  r  *  1,  ...»  m.  Then  X_(t)  and  X _ (t) 

"■  ""  ~  ■*—  —  n  n»T 

converge  in  distribution  to  I”ajXr (t)  and  Xr(t),  respectively. 

By  Corollary  4.7  we  conclude  that  Statement  (1.8)  holds. 

Theorem  4.8.  Let  tc(0,»)  and  r  *  1,  Then 


To  obtain  Statements  (1.9)  and  (1.10)  the  following  two  lemmas  are  needed. 


Lemma  4.9.  Let  k,  ae{l,2,..»}  and  a+k-1^  be  tbe 

order  statistic  of  a  sample  of  size  a  ♦  k  -  1  taken  from  the  population  U^. 

Then  ^^(a+q-l)”*!)  and  °a+jc_j  ^  are  identically  distributed. 

Proof.  The  spacing*  Vk-1,1'  °a+k-1.2  "  Vk-1,1*  ">  °a+k-l,k  '  Vk-l.k-1 

are  independent  exponential  rva’s  with  means  respectively  equal  to  (a+k-1)"*, 

(a+k-2)”*,  ...»  b~* ,  C Barlow -Proschan  (1975),  p.  59).  Hence  jjj  . (a*q-l)_1U 

H  1  H 

*n<1  ar*  '"'usl  ln  distribution.  <Vk-l,o'0)- 

Consequently  the  desired  result  follows.  || 

Lemma  4.10.  Let  t,  Be(0,«»),  and  r  *  1,  ...,  m.  Then  fim  E(Xn(t))P  = 

E{X(t))e  and  an  E(X„  „(t))P  »  E(X(t))8. 

n,r 


Proof.  By  Corollary  4.7  to  prove  the  results  of  the  theorem  it  suffices, 
(Chung (1974) ,  p.  95),  to  show  that  Sup  E{X  „(t)}P+1,  and  Sup  E{X  (t))^1  are 

®  ®  Sxl 

finite.  Since,  0£X  (t)£X  (t),  it  is  enough  to  show  that  Sup  E{X  (t)}p  <«. 

Il|4  *1  il 

Next  we  prove  the  preceding  statement.  Let  b  •  b^ ,  and  o  « 

Then  by  Lemmas  2.2  and  4.9  for  all  nc(l,2,...},  P(Xn(t)-b4k)  £  P{Jbalq“*U^£abt}  « 
■  (l-e"‘®bt)*t,  k  «*  0,  1,  ...  .  Further  for  all  ne{l,2,...),  E(Xn(t))P‘''1  • 

*  (B^l)/** (a+b)P  P(X  (t)-b>z)dz  £  (B+1)J*  ,0(<l*b+l)P(l-e”obt)<^  <  «*.  Consequently 
the  desired  results  follow.  || 


Clearly  from  Lemma  4.10  we  obtain  Statements  (1.9)  and  (1.10). 

Theorem  4.11.  Let  te(0,«*)»  and  r  »  1,  ...,  m.  Then  (i)  lim  EX _ (t)  ■ 

ITt  2a  t  a  t  n»r 

-  br(e  r  -1),  and  (ii)  tin  Var(Xn <r(t))  -  by(e  r  -e  r  ). 

Finally  from  Corollary  4.7  and  Letrnna  4.10  we  obtain  for  te(0,*>)  that 


lim  Var{Xn(t)}  = 


16 
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A  population  of  susceptible  individuals  exposed  to  m  contagious  diseases  is 
considered.  The  progress  of  this  epidemic  among  individuals  is  modeled  by  an 
m-dimensional  stochastic  progress.  The  components  of  this  process  represent  the 
number  of  infective  individuals  with  the  respective  diseases  at  time  t. 

A  class  of  m-dimensional  stochastic  processes  is  constructed.  These  processes 
describe  the  progress  of  the  epidemic  models  considered  in  the  sequel.  Exact  and 
approximate  formulas  for  the  Joint  and  marginal  state  probabilities  of  these  models 
are  obtained.  It  is  shown  that  the  approximate  formulas  are  vary  simple  functions  of 
time  while,  the  derivation  of  the  exact  formulas  Involve  tedious  computations. 
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